Flagellar dynamics of a connected chain of active, Brownian particles 
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Eukaryotic flagella are active structures with a complex architecture of microtubules, motor pro- 
teins and elastic links. They are capable of whiplike motions driven by motors sliding along filaments 
that are themselves constrained at an end. Here, we show that active, self-propelled particles that 
are connected together to form a single chain that is anchored at one end can produce the graceful 
beating motions of flagella. We use a combination of numerical simulations, scaling analysis and 
mean field continuum elastic theory to demarcate the phase diagram for this type of oscillatory 
motion as a function of the filament length, passive elasticity, propulsion force and longitudinal 
persistence of propulsion directions. Depending on the nature of the anchoring, we show that fila- 
ment either undergoes flagella-like beating or assumes a steadily rotating coiled conformation. Our 
system is simpler than its biological inspiration, and thus could be experimentally realized using a 
variety of self-propelled particles. 



Eukaryotic cilia and flagella are whiplike, elastic micro 
structures that undergo oscillatory beating to drive pro- 
cesses such as locomotion [1], mucus pumping [2], em- 
bryogenesis [3], and directed cell migration [4]. While 
the molecular mechanisms that control ciliary beating 
remain incompletely understood, it is well established 
that sliding forces generated by dynein motors attached 
to the microtubule-based backbone of cilia play a cru- 
cial role [5]. In addition to understanding how these 
active structures work in nature, there is growing inter- 
est in designing artificial analogs. Recent experiments 




FIG. 1: (a) (Clamped end) Left: Shapes taken by a stably os- 
cillating short filament over one beat cycle - the slight asym- 
metry is a consequence of noise. Right: A self-contacting, 
stably beating, long filament. (Pivoted end): Left: A stably 
rotating filament for small lengths similar to those observed 
in motility assays [16]. Right: A tightly coiled shape for large 
lengths, (b) Filament curvature as a function of time and the 
distance from the clamped end for (top) non-periodic beating, 
and (bottom) periodic buckling. Simulation movies are in the 
supplementary information [22] 
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on a minimal motor-microtubule system [6] demonstrate 
cilia-like beating, and artificial beating systems driven by 
external periodically varying electromagnetic fields have 
been synthesized [7]. However, internally driven struc- 
tures capable of controllable beating patterns have yet to 
be developed. In this letter we use simulations and the- 
ory to identify a different mechanism that results in con- 
trollable, internally driven flagella-like beating or steady 
rotation, in microstructures comprised of connected self 
propelled units. The resulting shapes and motions are a 
consequence of a fundamental instability of a filament un- 
dergoing tangentially-directed compressive forces - also 
termed follower forces. 

Our model filament is a two-dimensional chain of con- 
nected, self-propelled colloidal spheres. They may be 
experimentally realized, for example, as Janus beads in 
varying contexts [8-10] and can potentially be connected 
to form active filaments, either through face-face attrac- 
tions or by passive tethers [11-14]. Here we assume that 
the filament has N self-propelled polar spheres each of 
diameter a located at coordinates (i = l,..,iV). The 
spheres are connected by harmonic potentials of equi- 
librium length b, U\ = y (|r^ + i — — 6) 2 , and 
resistance to bending is implemented via a three-body 

bending potential Uh = | SiL^ 1 — b^j , where 

bi = (r^ — Ti-i)/\ri — Ti-i \ is the unit bond vector and n is 
the bending rigidity. Each sphere exerts a local force on 
the filament due to self-propulsion. This active force f* of 
constant magnitude acts on each sphere along its polarity 
vector p^, such that f* = f p pi with |p^| = 1. We consider 
self-linked Janus particles with a homogeneous polarity 
along the filament, i.e. the polarity vector itself tends 
to point along the linking (tangent) direction hi under 
the influence of a harmonic potential U a = (Pi — b^) 2 . 
For finite ft a , the angular potential biases the propulsion 
forces to act along the local tangent vector of the fila- 
ment; in the limit ft a — >• oc, the propulsion forces are 
perfectly aligned with the local tangent. The arc- length 
parameter s in the range < s < i = Nb parameterizes 
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the coarse-grained position along the filament. 

The dynamics of the filament is determined by the 
evolution of the sphere positions and orientations p^, 
which we simulate using overdamped Brownian Dynam- 
ics (i = 1, .., N), 
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Here F- = —dUi/dr- 1: = —dU^/dr^, and / p p are re- 
spectively the bond, bending, and propulsion forces, F a 
is the angular force derived from the angular potential 
?7 a , and Ff x = —dU ex /dri is an excluded- volume re- 
pulsive force given by a truncated and shifted Lennard- 
Jones potential for interparticle distance r: U ex = 

4e 
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e if r < 2e , and zero otherwise, and 

1/k^T. D and D T are translational and rotational 
diffusion constants, which in our low-Reynolds-number 
regime are set to satisfy the Stokes-Einstein relationship 
D r = 3D/a 2 . C and 

£ R are unit variance Gaussian white noise forces and 
torques respectively, each of which satisfies (Q) =0 and 
(Ci(t)Cj(t')} = 25 (t — t')8ij. As the viscous mobility of 
a sphere is isotropic, our simulations neglect differences 
in viscous resistance to normal and tangential modes of 
filament motion. However, simulations including full hy- 
drodynamic interactions, using a hybrid simulation tech- 
nique [15], display very similar qualitative features (see 
SI §111). We make the equations dimensionless by using 
<j and k^T as basic units of length and energy, and a 2 / D 
as the unit of time, and we set e = k^T. 

For our numerical simulations, runs were initialized us- 
ing a straight configuration with the filament vertically 
oriented along e x , with all p^ initially along — e x . One 
end of the filament corresponding to s = i was always 
free. The anchored end 5 = was either clamped ver- 
tically (Fig 1(a) top) [18] or attached to a frictionless 
pivot (Fig 1(a) bottom) with the filament free to rotate. 

Results for a clamped filament: For the clamped case 
b = e x , and r = 0. When the internal propulsion force 
/ p exceeds a critical value f c (£, « a ), the filament buck- 
les, in a manner similar to a self-loaded elastic filament 
subject to gravity [19]. However the post-buckled states 
are quite different (see SI-§I) since the direction of active 
force density follows polarity vectors p^ and thus tends 
to point along the filament axis. 

At fixed filament stiffness ft, and for f p > /c(A ft a )> 
we find that the magnitude of the polarization stiffness 
ft a controls the long time dynamics of the filament. In 
figure l(b - top) we show the local filament curvature as 
a function of time and arc length. For small values of 
ft a , thermal diffusion controls the local orientation of p^ 
and the propulsion activity is uncorrelated along the fil- 
ament. In this regime, the filament dynamics is marked 
by transients resulting from the bending generated by 
the particle propulsion, but no coherent patterns. When 
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FIG. 2: (a) Length of the end-end vector L ee for periodic 
beating (full line) and non-periodic buckling (broken line), (b) 
Power spectral density of the end-end distance for filament 
length £ = 80, propulsion force f p = 20, and angular force 
constant k& — 20, (c) critical propulsion force f c (£) for various 
filament lengths, at angular force constant K a = 100. The 
solid line is f c = Ck/£ 3 with C ~ 78, consistent with (2). (d) 
The beating frequency of filaments against propulsion force 
for filament lengths £ = 40 (•), £ = 50 (■), t = 80 (♦), and £ = 
100 (A). The dashed line corresponds to the scaling law(3). 
Parameters are dimensionless as described in the text. For 
very long lengths, the filament starts to self-contact (inset) 
resulting in a decrease in the exponent due to contact friction. 
The value of bending stiffness is k = 2 x 10 4 in all figures. 



K a is increased above a critical value the polarities of 
the spheres align strongly with the local tangent. The 
resultant self-propulsion force is strongly correlated with 
the filament tangent and the filament oscillates with pe- 
riodic, large amplitude wavy motions shown in Fig l(b - 
bottom) that propagate from the proximal (clamped) to 
the distal (free) end. This beating profile is very similar 
to that of flagella in eukaryotic cells though the under- 
lying physics differs fundamentally, since there is just a 
single filament made of active particles here. 

We quantify the regularity of these oscillations by 
m easuring the length of the end-end vector L ee = 
sj (r^y — ri).(r/v — ri), as a function of time, shown for 
two values of ft a , ft a < k\ and ft a ^> in Fig. 2a. While 
L e e displays large variations in time for both values of 
ft a because of the large propulsion force, the profile for 

^> ^ a is periodic, and its power spectral density (PSD) 
shows a distinct frequency maximum (Fig. 2b). 

To understand our numerical experiments and obtain 
estimates for the critical buckling load and the frequency 
of ensuing oscillations in the periodic regime, we consider 
the limit ft a — ^ oo, and coarse-grain the chain of spheres 
into a slender, elastic filament of length i and bending 
stiffness n. The force due to the self-propulsion trans- 
lates to a compressive force per unit length of strength 
/ p oriented anti-parallel to the local tangent vector. For 



the clamped filament, the straight filament becomes un- 
stable to time-dependent shapes beyond a critical value 
of / p , when the resultant internal propulsion force f p £ 
deflects the tip by a small transverse distance h <C £ 
leading to an effective filament curvature 0(h/£ 2 ). Bal- 
ancing moments about the base then yields hf p £ ~ Kh/£ 2 
and thence the critical force beyond which the straight 
filament is no longer stable, 



(a) 
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This scaling relationship closely agrees with computa- 
tionally determined values (Fig 2c). To determine the 
oscillation frequency we note that when the filament 
buckles, it is always under compression and no small 
amplitude, stable steady solution exists (SI-§1). When 
fp ^> /c? the characteristic length over which the active 
compression is accommodated scales as A ~ (ft/fp)^ <^ t~ 
In the over-damped limit, all the energy supplied by the 
self-propulsion transforms first into elastic bending en- 
ergy before being ultimately dissipated viscously. In a 
time cj _1 , the energy dissipated viscously is the prod- 
uct of the force per unit length t]±\uj, the characteristic 
deflection A, and the velocity uj\. This dissipation has 
to balance the active energy input into the system due 
to the self propulsion f p X 2 u and thus f]±X 3 u 2 ~ f p X 2 u. 
Using our previous expression for A ~ (^// P ) 5 , we obtain 

w-^x- 1 (3) 

which is also in excellent agreement with our simulations 
(Fig. 2d). A systematic derivation of the equations gov- 
erning the filament dynamics that builds on a local re- 
sistivity formulation [20] relating filament bending to its 
velocity yields two coupled non-linear equations for the 
tension in the filament T and (9, the angle the filament 
makes with e x (see SI- §2 for details) and confirms this 
simple scaling picture. 

Results for a filament anchored at a pivot. The beat- 
ing discussed thus far is a consequence of the clamped 
boundary condition which prevents rotation as well as 
translation. We next performed simulations for a fila- 
ment moving about a frictionless, pivoting end at s = 
such that i*o = whereas bo is unconstrained. 

For small values of ft a , and with the contour length £ 
and rigidity n held fixed, the filament end-end length L ee 
displays large irregular variations and the end-end vector 
L ee undergoes irregular rotation about the fixed end. In- 
creasing ft a results in the active forces being increasingly 
correlated along the contour. The post-buckled filament 
now assumes a steadily rotating bent shape and the value 
of L e e does not vary in time (Figure la). 

The rotation frequency extracted from simulations by 
calculating the orientation of L ee as a function of time is 
plotted in Fig 3 as a function of f p . For £ <C 37t(k/ f p )z , 
the wavelength of the bend is larger than the length 
of the filament and the steady rotational frequency ob- 
served varies with force in accord with equation (3). The 
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FIG. 3: (a) Frequency of steady rotation for a filament an- 
chored at pivot, with length £ = 160, angular force constant 
10 2 , and bending stiffness k — 2 x 10 4 . The in- 
set shows the radius of the coiled shape as function of f p . 
For very large forces, long filaments yield rotating closed coils 
and the frequency scales differently due to the excluded vol- 
ume constraint, (b) Transient shapes for very large f p as the 
filament coils and eventually undergoes steady rotation. 



shapes we obtain and the dependence of the frequency 
on the active force compare well with the experimen- 
tal observations of motility assays where filaments ani- 
mated by underlying molecular motors encounter pinning 
sites (defects) and start to rotate [16]. Beyond a critical 
value of / p (for £ fixed) or when £ ^> 37r(tt// p )3 (for 
/ p fixed), the wavelength of the bend becomes smaller 
than the length of the filament thus enabling self-contact. 
Since excluded volume interactions prevent overlap, the 
filament assumes increasingly coiled shapes as shown in 
Fig la and in Fig 3a. Once coiled, the frequency of steady 
rotation depends strongly on friction due to relative slid- 
ing and excluded volume forces. Our simulations indicate 
that for fixed £, the frequency almost scales as f 2 (Fig 
3a). Relaxing the excluded volume constraint results in 
completely overlapping coils and a frequency that scales 
as /p^ 3 , consistent with (3). 

Snapshots of the transient rotating shapes seen for 
long filaments (Fig 3b) show evidence of multiple time 
scales in the coiling process. Initially, the filament ro- 
tates and and adjusts its curvature until the free end al- 
most touches itself. At this point, the remaining length 
is accommodated in adjacently placed coils. Finally, as 
steady rotation is attained, the effective radius of the in- 
nermost coil R c reduces gradually, the coils tighten and 
the number of coils increases slightly. The radius R c de- 
creases very sharply with f p before eventually bending of 
the coil is balanced by the lateral forces due to excluded 
volume, resulting in a very slow decay with f p (Fig 3a 
inset). 

Dynamical correlations along filament: Having delin- 
eated the different dynamical regimes of the active fila- 
ment as a function of the active force and the anchoring 
condition, we now investigate the dynamical behavior of 
the connected self-propelled particles as a function of the 
angular stiffness ft a and self-propulsion force density f p . 
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FIG. 4: Phase diagram for (a) a clamped filament and (b) a 
rotating filament with bending stiffness k = 2 x 10 4 . Green 
triangles correspond to no beating, blue symbols to irregular 
beating, and red squares to regular beating. The blue dashed 
curve (valid for K a ^> 1) corresponds to equation (4). The con- 
stants were determined by eye and are (a) (A, B) = (85, 80) 
and (b) (A, B) = (36,70). 



In Fig 4, we show the phase diagrams for both clamped 
and pivoted filaments. 

For propulsion forces below the critical value f p < / c , 
we do not observe any statistically significant beating 
(green triangles) - the straight configuration remains sta- 
ble and variations in (L ee ) c± £ due to fluctuations in 
the directions of are insignificant. The critical value 
of f c that separates this regime from the unsteady (pe- 
riodic as well as non-periodic) regime is a function of 
the passive filament elasticity (£ and k) and also of 
Keeping ft a fixed and increasing f p to beyond the criti- 
cal value renders the straight configuration unstable and 
results in stable unsteady shapes. For small values of ft a 
(< 1) we observe irregular beating (Fig. lb, top). As 
ft a is increased (ft a > 1), enhanced correlation of par- 
ticle propulsion results in two trends - evident both for 
the clamped (Fig 4a) and pivoted (Fig 4b) cases. First, 
the critical force to destabilize the straight configuration 
decreases, saturating as k — )> oo to the continuum limit, 
consistent with (2). The increasing correlation between 
the self-propulsion directions of the spheres results in pe- 
riodic beating (red squares, Fig. lb, bottom). Interest- 
ingly however, at any ft a > 1, there persists a diminishing 
region of irregular beating between the no beating and 
the periodic beating regions. This behavior arises when 
the correlation time of propulsion directions is shorter 
than the filament oscillatory timescale in the perfectly 
aligned limit. 

To understand the location of the boundary between 
regular and irregular beating, we note that regular os- 
cillations will be interrupted when a critical number of 
neighboring spheres escape from an configuration driv- 
ing the beating to a configuration that is not unstable 
- chosen here to be the configuration with the polarity 
vector along the local tangent (and thus the force being 
extensional). The mean time required for this transition 
to occur T c may be simply estimated using a classical 
(mean) first-passage time analysis [21]. The angular po- 
tential in terms of the angle £ between p and the tangent 
vector is U a = ft a (l — cos£). The base of the potential 
corresponds to ( = and the energy when ( = tt is 2^ a 



(in k\)T units). For ft a > 1, equilibration in the angu- 
lar escape coordinate £ may be assumed with transitions 
occurring diffusively. Using the expression for the fre- 
quency of crossings in the over-damped limit (SI - §3), 
we find that the critical curve satisfies 



(4) 



Here A is a constant that ensures that the critical force 
matches the continuum limit (ft a — >> oo), and B is to be 
interpreted as the minimal size of the patch (equivalently 
force) required to cause beating irregularities. While 
strictly valid for ^> 1, this expression - plotted as 
the dashed blue line in Figs. 4a and 4b - satisfactorily 
captures the shape of the boundary for > 1. 

Conclusions. Our study suggests a simple proposal to 
mimic the beating of eukaryotic flagella. Rather than 
having motors that walk on adjacent filaments that are 
clamped at an end, we have shown that we can gener- 
ate actively oscillating and rotating filaments using con- 
nected self-propelled particles. The frequency of oscilla- 
tions and thereby the swimming speed and fluid forces 
can be controlled by varying the dimensionless parame- 
ters in our problem - the ratio of the chain stiffness to 
the polarity stiffness ft/ft a , the scaled active force £ 3 / p /k, 
and the aspect ratio of the chain £ /a, or equivalently 
the number of active particles. Furthermore, if the an- 
gular potential is written as U a = 4f[(p$ — b^) 2 — C], 
where < C < 1, thus introducing a preferred direction 
different from the local tangent, the initial symmetry- 
breaking bifurcation is eliminated and the filament will 
always beat. All of these parameters may be accessible 
experimentally. 

Recent efforts to manipulate connected passive colloids 
by electrical fields [11], microfluidics [12], lock-and-key 
type interactions [13] and heat [11, 14] have been success- 
ful in yielding externally actuated filaments with control- 
lable stiffness. Extension of such techniques using diffu- 
sophoretic Janus particles as templates should yield inter- 
nally controlled self-propelling filaments, just as a small 
variation of the motility assay for rotating filaments [16], 
by clamping an end, will lead to beating. In a broader 
context, in contrast to most studies of locomotion at low 
Reynolds number which prescribe the shape of the or- 
ganism (typically as a slender filament with prescribed 
kinematics), here we prescribe the active forces locally, 
and calculate the resulting shapes. If the anchored end 
does not have infinite resistance the self-propelled parti- 
cles will propel the whole chain, the study of which is a 
natural next step. 

After our work approached completion, we learned of a 
related numerical model [17]. However, our work differs 
in using a different formalism, different boundary condi- 
tions, as well as a consideration of both scaling relations, 
and the underlying mechanisms. 
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I. COMPARISON TO THE SHAPE OF A 
COLUMN BUCKLING DUE TO SELF WEIGHT 

In the continuum, noise-free limit one may treat the 
array of connected active swimmers as a thin, elastic, 
inextensible filament - valid provided Kb 2 jk^T ^> 1 and 
a jt <C 1. The filament bends due to the action of com- 
pressive, active forces and the constraint of inextensibil- 
ity. Such an active filament, clamped vertically at one 
end and free at the other, remains stably vertical for 
small force (compression) densities. At a critical value 
of the force density however, the straight shape is un- 
stable to lateral perturbations and yields to a buckled 
shape. 

Before buckling, the compressive forces all act verti- 
cally towards the clamped end and this scenario resem- 
bles the classical problem of a filament buckling under 
gravity. Originally proposed and solved by Euler in the 
18th century [1] this has been revisited again more re- 
cently by Keller and co-workers [2] (also see references 
therein). The correspondence arises due to the non- 
generic initial configuration with the filament clamped 
vertically. In the following, we will compare equations 
for this classical (passive) buckling problem with those 
for an (active) self-propelled filament and demonstrate 
that the post buckled states are very different. 

Consider a cylindrical, slender column of length i, di- 
ameter (j, cross-sectional area A = 7rt 2 /4, bending stiff- 
ness n and mass density p clamped vertically upwards to 
a rigid flat surface at its base 5 = and free at the other 
end s = £. No external forces or torques act at the free 
end. We choose the origin of our Cartesian coordinate 
system to coincide with 5 = 0, the x axis to point verti- 
cally upwards and gravity to point downward such that 
g = —ge x . We parameterize the shape of the column by 
the angle, 6(s,t), the centerline makes with e x and de- 
compose the force resultant at a point s along its length 
in terms of its Cartesian components, F = F x e x + F y e y . 
Torque and force balances on a elemental length yield 
(primes henceforth denoting differentiation with respect 
to arc-length, s) M" = F x 6' cos6+F y 6' sinO+Apg sin (9, 
F' x = Apg, and F' y = 0. Since there are no forces acting 
along the y direction, we may integrate to find F y = 
and F x = Apg(s — £). Substituting this in the equation 
for the moment and using the definition M = k6' \ we find 
that the filament shape satisfies the differential equation 
k6'" = Apg((s — £) sin#y. In the small deformation limit 
sin# & 6, we scale length with i and we obtain an equa- 



tion for the angle 

0"' + (a (1 - s) 0)' = 0, a= (1) 

with the associated boundary conditions 6(0) = 0, 
<9'(1) = and 6" (I) =0. 

What is the corresponding equation for the filament 
of connected self propelled particles? Assuming that a 
static steady solution exists, we use a moment balance 
as before. While it is easier to work in a local refer- 
ence frame that uses the tangent and normal vectors 
at any point as the basis vectors, we choose to retain 
the Cartesian nomenclature to allow for comparison with 
equations for buckling due to gravity. The torque and 
force balances now yield k6"' — 6'(F X cos 6 + F y s'm6), 
F' x = f p cos 6 and F y = f p sm6. Now making the small 
deformation assumption and using the fact that no exter- 
nal forces are present at the free end, we find that these 
equations reduce to 

6"' + [3(1 - s)6' = 0, p = ^- (2) 

where the constant f3 = (f p £ 3 /^) - a form different from 
equation (1). 

The next question to answer is if steady small defor- 
mation solutions exist in either of these cases. We seek 
non-trivial solutions in terms of the Airy functions A 
and B. Let us first consider the column buckling due to 
self-weight for which a solution %i = A(X) + ^= B(X) 
exists with constant C undetermined at linear order and 
X ee as (s — 1). We use the expansions of the Airy func- 
tions in terms of Bessel functions J, i and J_i for s > 

^3 3 

[3] and determine the particular solution satisfying the 
boundary conditions and find that the critical buckling 

load a c is given by the implicit equation J_i ^|a^ =0. 

Thus there exists a non-trivial small deformation solu- 
tion. 

For the active filament with boundary conditions 
6(0) = 0, 6' (I) = and 6" (I) = 0, we again seek solutions 
in terms of the Airy functions. Recasting the equation 
for 6 as 6 = J Z(s') ds' , and defining (p = ft* (s — 1), we 
obtain the differential equation 

Z w -<pZ = 0, Z(0)=0, Z (p (0) = 0. (3) 

Solutions of the form Z — C\A(tp) + Ci1&(}p) can satisfy 
the required boundary conditions only when C\A(0) + 
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FIG. 1: Schematic of the bent state and reference sketch 
showing the force and torques acting on a elemental slice. 



C 2 B(0) = 0, and dA'(0) + C 2 B'(0) = are satisfied 
simultaneously. This is possible only when C\ — and 
C 2 = 0. One can have small deformations solutions that 
are unsteady or large deformation steady solutions (the 
latter being ruled out by our simulations) . 



II. COARSE GRAINED CONTINUUM MODEL 

FOR k~ 1 -> 

While the previous analysis suggests that small defor- 
mation static shapes are unstable, unsteady (oscillatory) 
or unstable solutions can still exist. To derive the equa- 
tions governing unsteady shapes that may arise, we be- 
gin by choosing a convenient parametrization for the fil- 
ament shape. Figure 1 is a schematic of the geometry 
and also shows the free-body diagram of the forces and 
torques acting on an elemental length of the filament. 
We coarse-grain the discrete filament of attached self- 
propelled Brownian particles into a continuous inexten- 
sible filament of length i and diameter a moving in the 
x—y plane. The active compressive forces are also coarse- 
grained into a force density f p acting anti-parallel to the 
tangent vector. Choosing the arc- length s as our variable, 
we locate the origin of our stationary co-ordinate system 
at the clamped base s = 0. The other end s = t is free to 
move in the x — y plane. The beating is characterized by 
the sequence of shapes generated as a test material point 
at s with cartesian coordinates (x p (s),y p (s)) moves in 
the Newtonian liquid of viscosity fi. In the long aspect 
ratio limit tja ^> 1, the angle 0(s,t) made by the center- 
line of the inextensible filament with the x axis serves as 
a convenient indicator of the filament shape. With this 
parametrization, we can express the dual vectors tangent, 
t(s, t) and normal n(s, t) to the centerline solely in terms 
of 0(s,t). Note that at (x pi y p ), the increments along the 
cartesian directions are related to ds by, dx = cos ds 
and dy = sin 6 ds. 

We now work in the small deformation limit sin# « 



and only consider forces that generate curvatures that 
are much greater than a. Referring to the free body dia- 
gram in Figure (1), we note that in the physically relevant 
non-inertial (small Reynolds number) limit, the viscous 
forces, elastic and active forces are all in balance. To 
make progress we decompose the total force resultant at 
a cross-section s, F(s, t) into its tangential T and normal 
TV components 

F = Tt + Nn. (4) 

Furthermore, we assume that the viscous forces are given 
by the values obtained from local resistivity theory [4] - 
thus the viscous force per unit length at s is 

f v = — (77||^||t + r]±u±n) (5) 

where m is the effective viscous resistance per unit length 
for motion of the filament along the tangent, rj± is the 
resistance per unit length for motion along the local nor- 
mal, un is local velocity of the centerline along the tan- 
gent vector and u± is the filament velocity along the nor- 
mal. The force balance then is 

- F(s) + F(s + ds) - f p t ds = -f v . (6) 

Using 

F' = (T' - NO') t + (N' + TO') n 

we balance force components to obtain the equations 
(TV' + TO') = rj ± u ± and (T' - N0 f ) = rj in + f p . A bal- 
ance of moments acting on the differential element yields 

M' + tx (Nn + Tt) = (7) 

Combining all the above expressions, we can eliminate 
TV and relate the bending moment per unit length to the 
angle using M = n0' to finally obtain the coupled non- 
linear equations for the tension, T and angle 

- k0"' + TO' = f] ± u ± (8) 

T' + K0"0'-f p =v\n (9) 

To close these equations, we need to relate the veloc- 
ity of the filament to its shape and properties. For an 
inextensible filament, we have following the steps in [4] 

ds „ , dt dO 80 
— = 0, and — = — n = — n 
dt dt dt ot 

and at the same time, 

f = (u\ l -e'u ± )t + (u' ± + e'u ll )n. 

Using these expressions to eliminate the filament veloci- 
ties in favor of the filament shape and shape changes and 
indicating derivatives with respect to time as subscripts, 
we find 

t = u' ± + u\\0', and = u\\ - u±0' (10) 
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and consequently 

T" = -«(0"0')' + (f P Y + ?l 0' (-K0'" + TO') (11) 

>/- 

and 

77_L0 t = -/c0"" + (re 7 ) 7 + — 0' (T' + ^0"0' - /p). (12) 

^11 

A full numerical solution to these equations under the 
constraint r]± = rj\\ will yield the pre-factors in equations 
(2) and (3) of the main text. 

III. A FIRST PASSAGE TIME CALCULATION 
FOR FINITE n a 

We adopt a simplified picture of the the local dynam- 
ics and fluctuating forces generated by the self-propelled 
spheres. Since our initial condition is a vertical filament, 
this special orientation results in a critical load for buck- 
ling before any sustained bending can occur. 

Let the total active force f p (s) be separated into a 
quasi- static average part and a rapidly fluctuating part: 
fp( s ) = (fp( s )) + fp( s )- Fluctuations arise as the ori- 
entations of individual active elements escape from a fa- 
vorable configuration (anti-parallel to the tangent) to an 
unfavorable configuration where the active force contri- 
bution is negligible - chosen to be along the normal, n(s). 

In the diffusion (noise) dominated limit K a « 1, the ac- 
tive, polar spheres are very animated, re-orienting them- 
selves constantly with respect to their neighbors and over 
very small time scales compared to the response of the fil- 
ament. Thus, we then have both K a «l and the charac- 
teristic filament beating timescale Tf ~ nf p 3 ^> D^ 1 . 
Irrespective of the magnitude of the constant part of f p 
the filament does not undergo correlated beating in this 
limit. In agreement with this statement, the simulations 
show only irregular beating or no beating for n a < 1 
(Figure 4 of the main text). 

Consider now the case when 1 <C ft a < oc. Simu- 
lations (Figure 4 of main paper) demonstrate three dis- 
tinct regimes. For low propulsion forces below the critical 
value / p < / c , we do not observe any statistically sig- 
nificant beating (green triangles in Figure 4 of the main 
text) - the straight configuration remains stable and vari- 
ations in (L e e) — ^ due to fluctuations in the directions 
of pi are insignificant. As f p increases at constant n ai 
we go from the no-beating regime (A) to an intermediate 
irregular beating regime (B) and finally to the regular 
beating regime (C). Note that the critical active force 
density that separates both transitions (A — > B as well 
as B — » C) depends on the passive filament elasticity (£ 
and k) and also on n a . 

To understand the physical mechanisms underlying 
these transitions, we first note that the scaling for f c 
derived in the main text (Eqn. 2) corresponds to the 
noise- less regime = 0. In the weak noise limit, the 



effective time to escape is extremely large due to the high 
barrier to rotation of polarity vectors. Thus for f p < / c , 
the filament is straight and for f p > f p the filament first 
buckles and then beats periodically with a frequency that 
scales as f p ^ 3 . 

With this in mind, consider fixing the active density f p 
to a value greater than f c and varying K a (moving along 
the horizontal axis in the figures). As n a is increased 
from unity, the enhanced correlation of particle propul- 
sion directions results in two trends - evident both for the 
clamped (Fig 4a) and pivoted (Fig 4b) cases. First, at 
some value of n ai the force density becomes large enough 
to destabilize the straight configuration (regime A) re- 
sulting in a regime (B) characterized by erratic beating. 
This /%-dependent critical density decreases as n a in- 
creases, saturating to the noiseless continuum value as 
k — » oo. In regime (B), the mean active force density 
is large enough to cause buckling over the entire length 
of the filament £ and the orientation potential well has a 
depth (proportional to K, a ) which is large enough that few 
polar spheres reorient. Thus, with increasing n a there is 
increasing correlation between the sphere self- propulsion 
directions as well as an increasing time over which the 
correlations are sustained. Since the self-propulsion di- 
rections tend to point along the filament tangent, there 
is a typical time, T c , for which the propulsion direction 
remains along the local tangent vector. The last trend is 
particularly relevant to whether or not there is sustained, 
regular, oscillatory beating. In particular, we hypothe- 
size that irregular beating (regime B) arises whenever 
the correlation time T c of propulsion directions is shorter 
than the filament oscillatory timescale in the noiseless, 
perfectly aligned limit, i.e. Tf ~ T]±( K /fp) 5 • 

In summary, irregular beating arises when T c <C Tb 
while regular beating (regime C) arises when Tb ^> T c . 
Thus we expect the critical curve characterizing the B — >> 
C transition to satisfy Tf ~ T c . 

The final part of the analysis requires estimating this 
correlation time T c . We first rewrite the harmonic angu- 
lar potential, U a = ^(pi — b^) 2 (which biases the polar 
spheres to orient anti-parallel to the tangent), in terms 
of the angle (. Thus, U a = ft a (l — cos^). We recognize 
that T c also corresponds to the mean passage time for 
a critical number of self-propulsion directions within a 
local patch of spheres to escape from a perfectly aligned 
configuration, in which they drive filament compression, 
to a stable non-compressive configuration. For the calcu- 
lation we choose the latter to be a configuration in which 
the polarity vector and the local tangent are in the same 
direction (£ = tt) so that the forces are extensional and 
thus no compressive instabilities are possible. 

The energy difference between the two configurations 
( — and ( = 7r respectively is 2^ a . For n a > 1, equi- 
libration in may be assumed, with the polar spheres 
crossing the barrier diffusively. In the over-damped limit, 
this flux can be related to a frequency of crossover or the 
disrupting frequency. While a complete calculation of the 
escape frequency is complicated, an asymptotically accu- 
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Green = w/ HI 
Red = w/o HI 



FIG. 2: (a) Trajectories of material points on the filament 
with a clamped boundary condition, with (blue) and without 
(red) full hydrodynamic interactions (HI). Inset: the end-to- 
end distance L ee of beating filaments as a function of time, 
(b) The configurations of a filament with a pivoting boundary 
condition for simulations with (green) and without (red) full 
HI. For (a) and (b) the filament length is £ = 80, and the 
active force density is f p = 10. (c) Schematic of a filament 
configuration where relative motions of adjacent filaments en- 
counter increased viscous drag when full HI are considered. 



rate estimate obtained using the steepest descent method 
and valid in the over-damped, viscosity controlled limit 
gives [5]: 

r ~ ^ e - 2 ^ 

27T7b 

where 7b depends on the resistivity in the vicinity of the 
escape point (£ = 7r), the frequency of the stable basin 
is given by Uq ~ |£7"(0)|2 ~ and the frequency of 
the barrier is given by ~ \U"(tt)\^ ~ Since 
the diffusivity and thus the viscosity are constant in all 
our simulations, we ignore multiplicative factors. Setting 
this frequency to be of the same order as the zero-noise 
frequency (Equation 3 in the main text) we obtain a pre- 

- 9 3 

liminary estimate for the critical line of / p ~ k£ (e~ a ) 4 . 

This estimate needs to be modified as follows. First, 
a non-zero critical density for buckling arises due to the 
non-generic initial configuration. However once buckled, 
the net force determines the beating and the beat fre- 
quency is a function of / p , not f p — f c . Thus we in- 
troduce A*, an adjustable parameter which ensures that 



the critical force tends to the noiseless continuum limit as 
ft a — » oo. Second, the irregular beating behavior requires 
that a critical number of motors escape from the well, 
which we account for by introducing a pre-factor B* as a 
free parameter. Satisfying these conditions leads to the 
estimate 

f p ~A*+B* K I(e- 2K *)i. (13) 

We rewrite this in dimensionless form by redefining A = 
A*£ 3 /k and B = £ 3 / n to obtain equation 4 of the main 
text. 

^~A + BKl(e- 2 ^)i, (14) 



IV. EFFECTS OF HYDRODYNAMIC 
INTERACTIONS 

In simulations that were discussed in the main text, 
we used a freely draining model for polymer, where 
anisotropy of drag on the filament as well as hydrody- 
namic interactions between different parts of the filament 
due to the flow set up in the surrounding fluid were ne- 
glected. 

In order to capture the effect of full hydrodynamic in- 
teractions (HI) on the dynamics of attached filaments in 
the limit ft a — >• oo, we used a hybrid simulation tech- 
nique, in which molecular dynamics simulations for the 
filament were combined with a mesoscale hydrodynamic 
simulation method called multi-particle collision dynam- 
ics (MPC) for solvent [6]. In this approach we model the 
solvent as a collection of TV point-like particles of mass m, 
whose velocities are determined by a stochastic process. 
Two steps are performed at each time point to evolve a 
trajectory. In the streaming step, the particles move bal- 
listically for a time interval h that may be understood as 
a mean collision time. In the second step (the collision 
step), the particles are sorted to cells of a square lattice 
with lattice constant a, and the particle velocities rela- 
tive to the center-of-mass velocities of the cell, are rotated 
by an angle a. The direction of rotation is chosen ran- 
domly. The conservation of mass and momentum in each 
cell yields accurate long range HI. The dynamics of the 
active filament is meanwhile simulated using a standard 
velo city- Ver let algorithm, and the filament-fluid interac- 
tion is implemented by including the filament monomers 
in the collision step and allowing for the appropriate mo- 
mentum transfer. 

Simul ation para meters for the solvent were a = 130°, 
h = ^/ma 2 /^T, and the mean number of particles 
per cell (TV) = 20. For the filament, the monomer 
mass M = m(7V), the bond length was b = a, and the 
monomer diameter was a = b. A constant temperature 
was maintained by local rescaling of the solvent velocity. 

We next compared the dynamical behavior of both 
clamped and pivoting filaments with and without full HI. 
For the clamped boundary condition, we found the same 
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critical active force density for buckling f c as for the sim- 
plified local hydrodynamics model, and for f p > f c the 
filament interacting via full HI displays the same flagella- 
like beating as does the simplified model. Furthermore 
the expression for beating frequency as a function of ac- 
tive force density f p is the same in both cases. A detailed 
comparison of the beating patterns of the two predictions 
is made by tracking specific material points along the 
filament contour. We plot in Figure 2(a) the trajecto- 
ries of the three material points located at s = £/4, £/2 
and i when i — 80. The results demonstrate that beat- 
ing patterns in the presence and absence of full HI are 
qualitatively similar, with the model filament transcrib- 
ing a figure-of-eight in both cases. We see that including 
full HI leads to slightly smaller lateral amplitudes (Fig. 
2(a)) and increases the beating frequency (inset). The re- 
duction in both relative motions between filament parts 
as well as a reduction in beating amplitude (due to the 
increased viscous interactions) when combined with the 
constant energy input due to activity may account for 
the increase in frequency. However, we reiterate that the 
differences in beating patterns with and without HI are 
not qualitative and that expressions for the frequency in 
terms of the active force density, f p and the critical force 
for buckling are the same. A pivoting boundary condi- 
tion with full HI yields a rotating filament with similar 
shapes as found for the simplified model (Figure 2(b)). 



Non-local hydrodynamics play a more important role 
for very long filaments, where parts of the filament sepa- 
rated by distances much greater than the filament thick- 
ness a come close to each other (Figure 2c) . Our minimal 
simplified model incorporates self- avoidance accurately 
via the excluded volume interaction force, and isotropic 
hydrodynamic resistivity tensors to compute the local 
viscous drag. A slightly more accurate expression for the 
drag per unit length f v (s) at material point x with veloc- 
ity xt that incorporates the effect of anisotropy (thus the 
filament is no longer a freely draining polymeric filament) 
is given by 

f„(M) = -^(2^-tt).x, 

This expression however still neglects the full non-local 
hydrodynamics that arises naturally from our hybrid 
MPC simulation. While the qualitative response (beat- 
ing or rotating) is expected to be the same for long fila- 
ments, the increased role of viscous forces when the fila- 
ments are very close to one another (as in highly looped, 
beating configurations or tightly wound rotating shapes) 
modifies the beating patterns significantly. A systematic 
analysis of the changes in beating frequencies and shapes 
of long, self interacting filaments, due to fluid mediated 
interactions is a topic for our future study. 
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